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Abstract. We describe a new code for the deep anal- 
ysis of stellar fields, designed for Adaptive Optics (AO) 
Nyquist-sampled images with high and low Strehl ratio. 
The Point Spread Function (PSF) is extracted directly 
from the image frame, to take into account the actual 
structure of the instrumental response and the atmo- 
spheric effects. The code is written in IDL language and 
organized in the form of a self-contained widget-based ap- 
plication, provided with a series of tools for data visualiza- 
tion and analysis. A description of the method and some 
applications to AO data are presented. 

Key words: instrumentation: adaptive optics - method: 
data analysis - techniques: image processing - techniques: 
photometric - astrometry - stars: imaging - galaxy: center 



1. Introduction 

StarFinder is a code developed within the frame of the 



ESO PAPAO (Currie et al. |2000aj gOOOb|) program for 
the reduction of Adaptive Optics data. It was originally 
designed (Diolaiti et al. 1999) to analyze AO images of 



very crowded fields, like the PUEO frame of the Galactic 
Center shown in this paper, which contains about 1000 
detectable stars in a field of view of 13x13 arcsec^. This 
image is an example of an AO observation: the PSF is 
Nyquist-sampled and characterized by a complex shape, 
showing a sharp peak, one or more fragmented diffraction 
rings and an extended irregular halo. Moreover, due to 
the small field of view, the imaged region is approximately 
isoplanatic and the PSF may be considered space invari- 
ant. Under the assumptions of isoplanatism and Nyquist- 
sampling, StarFinder models the observed stellar field as 
a superposition of shifted scaled replicas of the PSF lying 
on a smooth background due to faint undetected stars, 
possible faint diffuse objects and noise. 

The procedure derives first a PSF digital template 
from the brightest isolated field stars; then a catalogue 



of presumed objects is formed, searching for the relative 
intensity maxima in the image frame. In the following step 
the images of the suspected stars are analyzed in order of 
decreasing luminosity. In this phase a catalogue includ- 
ing the accepted objects is formed and a synthetic image 
of the observed field is constructed, placing an intensity- 
scaled PSF template in the position of each identified star. 
Each suspected object in the original list is accepted on the 
basis of its correlation coefficient with the PSF template; 
the relative astrometry and photometry of the source are 
determined by means of a fit, taking into account the con- 
tribution of the local non-uniform background and of the 
already detected stars, known from the synthetic field. As 
the analysis proceeds, fainter and fainter sources can be 
successfully analyzed, discriminating their central peaks 
from the secondary bumps of light in the distorted diffrac- 
tion rings of the neighboring more luminous already iden- 
tified stars: in this way the synthetic field becomes more 
and more similar to the observed image. Residual unex- 
plained features of the image may be further analyzed 
checking for indications of blended groups, at separations 
smaller than the PSF FWHM. 

StarFinder should be regarded as a tool for high pre- 
cision astrometry and photometry of crowded fields ac- 
quired under the above assumptions of accurate PSF 
knowledge, isoplanatism and correct sampling. An appli- 
cation of StarFinder to HST undersampl ed da ta handled 
by dithering techniques (Fruchter ct al. 1997) has been 
shown in Aloisi et al. ( 200C| ). In this paper it can be seen 
that the results obtained by our method are comparable 
to those obtained by DAOPHOT (Stetson |1987|). 



As far as the PSF is concerned, the accurate knowledge 
of its features outside the central peak is fundamental to 
perform a deep study of a stellar field, achieving accu- 
rate photometry of faint stars and avoiding either false 
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detections or star loss (Esslinger & Edmunds 1998). If the 
PSF template cannot be extracted directly from the field, 
due to extreme crowding or lack of bright isolated stars, 
StarFinder can still be applied using a PSF estimated by 
means of other methods, as the reconstruction technique 
proposed by Veran et al. (1997) for AO observations or the 
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TINY TIM simulation software for HST (Krist & Hook 



1999) 



Much more intriguing and difficult to solve is the case 
of a field with space variant PSF, either due to anisopla- 
natic effects as in AO observations or to design and control 
in HST images. In general the analysis of an anisoplanatic 
field requires the knowledge of the local PSF. A method 
to reconstruct the off-axi s PSF in AO imaging has been 
proposed by Fusco et al. ( 2000 ). StarFinder, in its present 
version, can analyze frames with space invariant PSF or 
sub- frames in which the isoplanicity condition is nearly 
satisfyied, as will be shown in Section 4; the extension 
to the space variant case is in progress, with preliminary 
results presented by Diolaiti et al. ( 2000a , 2000b ). 

The paper is organized as follows: the general features 
of the algorithm are described in Sect. 2; Sect. 3 deals 
with more technical aspects and might be skipped on a 
first reading; the method is validated on simulated and 
experimental data in Sect. 4; details on the IDL code are 
presented in Sect. 5; Sect. 6 includes our conclusions and 
future plans. 
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Fig. 1. Left: sub-image extracted from a simulated stellar 
field. The crosses indicate the objects within the region 
of interest: a central star (to be analyzed in the current 
step), a brighter source (already known), a fainter one (to 
be examined later) and the PSF feature of a much brighter 
star, represented by the structure in the upper-left part of 
the sub- image. Right: corresponding sub-region extracted 
from the stellar field model, containing one replica of the 
PSF for each star detected so far. 



2. Code description 

2.1. PSF estimation 

The accuracy of the PSF estimate is primary in 
StarFinder, since the PSF array is used as a template for 
all the stars in the image to be analyzed. The user selects 
the most suitable stars, which are background subtracted, 
cleaned from the most contaminating sources around, cen- 
tered with sub-pixel accuracy, normalized and superposed 
with a median operation. The centering is performed by 
an iterative shift of the stellar image in order to cancel the 
sub-p ixel offset of its centroid (see Christou & Bonaccini 
19961 ). The median operation, which is performed pixel- 



by-pixel, is preferred to the mean because it is less sen- 
sitive to anomalous pixels (outliers) which might appear 
in one or more stellar images among the selected ones. 
The retrieved PSF is post-processed in order to reject any 
residual spurious feature and to smooth the noise in the 
extended halo. 

The PSF extraction procedure also reconstructs the 
core of saturated stars by replacing the corrupted pixels 
with the central part of the PSF estimate. Accurate po- 
sitioning is achieved by means of a cross-correlation tech- 
nique, while the scaling factor is determined with a least 
squares fit to the wings of the star to repair. For a detailed 
description of the procedure, see Sect. 3.5. 

2. 2. Standard analysis of a stellar field 

At first we build a list of objects, the presumed stars, 
which satisfy the condition 



i{x,y) > b{x,y) +t 



(1) 



where i{x,y) is the observed intensity, b(x,y) the back- 
ground emission and t a fair detection threshold. The pre- 
sumed stars are analyzed one by one by decreasing in- 
tensity. To illustrate a generic step of the analysis, we 
consider the (n + 1) — th object in the list, after the ex- 
amination of the first n. A small sub-image of fixed size is 
extracted around the object (Fig. |^, left). This sub- image 
may contain brighter stars formerly analyzed, fainter ob- 
jects neglected in the current step and features of other 
stars lying outside the sub-image. The information on the 
brighter sources is recorded in a synthetic stellar field 
(Fig. ^ right), defined as the sum of two terms: a su- 
perposition of PSF replicas, one for each star detected up 
to this point, and an estimate of the background, assumed 
to be non uniform in general. The local contribution due 
to the brighter stars and the background, derived from 
the synthetic field, is subtracted from the sub-image. If 
a statistically significant residual remains, it is compared 
to the central core of the PSF by means of a correlation 
check. If the correlation coefficient is higher than a pre- 
fixed threshold then the object of interest is rated similar 
to the PSF and accepted. The accurate determination of 
its position and relative flux is attained by means of a local 
fit, in which the observed sub- image is approximated with 
the multi-component model described in Sect. 3.6. The ac- 
tual size of the fitting region is comparable to the diameter 
of the first diffraction ring of the PSF. This choice ensures 
that the information represented by the shape of a high- 
Strehl PSF is included in the fitting process to achieve 
better accuracy and prevents the growth of the number of 
sources to be fitted together. For the central object of our 
example a single component fit is performed and the con- 
tribution due to the brighter stars is considered as a fixed 
additive term. A multi-component fit is performed when 
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Fig. 2. Flow-chart of the algorithm for stars detection and 
analysis. 



the star is in a very compact group, at separations com- 
parable to the PSF FWHM. If the fit is acceptable, the 
parameters of the new detected star are saved and those 
of the already known sources, which have been possibly 
re-fitted, are upgraded. The new star and an upgrade of 
the re-fitted sources are added to the synthetic field. 

This analysis is performed for each object in the initial 
list (a flow-chart illustrating the operations of StarFinder 
is in Figure g) . To achieve a better astrometric and pho- 
tometric accuracy, at the end of the examination of all 
the objects, the detected stars are fitted again, this time 
considering all the known sources. This step may be iter- 
ated a pre-fixed number of times or until a convergence 
condition is met. 

At the end of the analysis, it is possible to stop the 
algorithm or instead perform a new search for lost ob- 
jects removing the detected stars and using an upgraded 
background estimate. It should be stressed that this image 
subtraction is just a tool to highlight significant residuals. 
Any further analysis is performed on the original frame, 
in order to take into account the effects arising from the 
superposition of the PSFs of neighboring sources. Gener- 
ally, after 2-3 iterations of the main loop, the number of 
detected stars approaches a stable value. 

2.3. Crowding and blending effects 

A binary star with different separation values (Fig. g) has 
been simulated to show how the code works with crowded 
sources. With a separation of 2 PSF FWHM the two com- 
ponents are well separated and the code analyzes them 
with the standard procedure described in the previous 
sub-section. In the other cases (separation from 1 to 0.5 
PSF FWHM) the secondary component is not detected as 
a separate relative maximum and it is lost. However, if the 
separation is not as small as 0.5 FWHM, a further itera- 
tion of the main loop enables the algorithm to detect the 
fainter component by subtracting the brighter one. This 
strategy forces the two stars to pass the correlation test, 
the principal component as a presumed single object and 
the secondary in a subsequent iteration of the loop. In a 
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Fig. 3. Simulated binary stars at various separations. 
From left to right, top to bottom the separation is 2, 1, 
0.75, 0.5 times the PSF FWHM. For aU the images the 
flux ratio is 2:1. 

way the iteration of the main loop is a de-blending strat- 
egy, because it enables the algorithm to detect stars whose 
intensity peak is not directly visible in the observed data. 
This strategy fails when: 

— the binary is very noisy and the two components have 
similar flux. In a similar situation the principal compo- 
nent may not pass the correlation test, preventing also 
the detection of the secondary in a further iteration. 
These noisy blends may be recovered at the end of the 
overall analysis. 

— The two components are almost equally bright and 
have a separation close to the lower limit (1/2 PSF 
FWHM) . The residual corresponding to the secondary 
after subtracting the principal component may have 
a distorted shape and might not pass the correlation 
check. 

The latter case may be handled by a method based 
on a thresholding technique. The object is cut at a pre- 
fixed level, about 20% below the central peak, and trans- 
formed to a binary array, setting to 1 all the pixel above 
the threshold and to the pixels below. If the area of 
the pixels with value 1 is more extended than the PSF, 
the object is considered a blend and the secondary star 
may be detected by subtracting the brighter one; then a 
two-components fit allows accurate astrometry and pho- 
tometry of the two sources. This strategy can be itera- 
tively applied to multiple blends. It should be stressed 
that the area measurement is not reliable when the value 
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of the cutting threshold is comparable to the noise level: 
for this reason the de-blending procedure is applied only 
to objects with a suitable signal-to-noise ratio. Moreover 
the area measurement is reliable when the data are ade- 
quately sampled. This de-blending procedure is applied at 
the end of the last iteration of the main loop, when all the 
resolved sources have been detected: in this way the prob- 
ability that a single object may appear artificially blurred 
because of the contamination of still unknown sources is 
largely reduced. 

In a normal case, like the simulated field of Sect. 4.1, 
two or three iterations of the main loop find almost all 
(~99%) the stars that StarFinder may detect. The de- 
blending procedure described above adds ~1% more stars, 
without additional false detection. The number of lost ob- 
jects belonging to the first category described above is 
negligible (<<1%). 

Normally we perform two or three iterations of the 
main loop and apply the de-blending strategy only in very 
crowded fields. 



3. Technical aspects 

3. 1 . Bad pixels 

StarFinder includes an automatic procedure to repair 
known bad pixels, which are replaced with the median of 
the good data points in a suitable neighborhood. The use 
of this procedure is particularly recommended in the ex- 
traction of the PSF. Indeed the stellar images, which will 
be combined to form the PSF, are centered with sub-pixel 
accuracy by means of an interpolation technique and the 
presence of bad pixels may produce unpredictable interpo- 
lation errors. Even if replaced, the bad pixels are excluded 
from any further computation concerning the correlation 
coefficient and the fitting process. 

3.2. Background estimation 

A reliable estimate of the background is necessary to de- 
fine the detection condition and to compute the correla- 
tion coefficient of the presumed stars with the PSF. 

A straightforward technique to estimate the back- 
ground is to smooth the image by median filtering, re- 
placing each pixel with the median computed over a suit- 
able neighborhood, of size comparable to the characteris- 
tic width of a stellar image. This method tends to over- 
estimate the background underneath strong peaks. A more 
accurate approximation, even in crowded fields, is de- 
scribed in Bertin & Arnouts ( 1996| ). The image is par- 
titioned in sub-regions arranged in a regular grid and a 
local estimate is calculated for each sub-patch by means 
of an IDL im plem entation of the DAO PHO T SKY rou- 
tine (Stetson 1987), due to Landsman ( |1995D . This array 
of sky measurements is smoothed by median filtering and 
interpolated onto the same grid of the input image. 



It should be stressed that the background computation 
is unavoidably affected by the presence of bright sources. 
In general a more accurate estimate can be obtained after 
the analysis of the stellar field, when the most contami- 
nating sources are known and can be subtracted. 

3.3. Noise estimation 

The estimate of the noise is useful to define the detec- 
tion threshold and to compute the formal errors on the 
retrieved astrometry and photometry. 

The overall effect of photon and instrumental noise can 
be computed if the required parameters (detector gain, 
read-out-noise, dark current, etc.) are known. Otherwise 
an estimate of the mean background noise can still be 
obtained by means of histogram fitting techniques (Al- 
moznino et al. 1993, Bijaoui 198C| ). Assuming that the in- 
tensity of the sky radiation is distributed normally around 
a typical value, the histogram of the observed intensity 
levels should be quite similar to a gaussian distribution, 
whose mode and standard deviation represent respectively 
the sky value and the associated noise. Actually the con- 
tamination due to stellar sources produces a high-intensity 
tail and an artificial broadening of the histogram, which 
prevents an accurate estimate of the background noise. 
This problem can be partially overcome by removing the 
signal from the image, leaving only the pixel-to-pixel vari- 
ations associated to pure noise: a reasonable estimate of 
the signal to subtract for this purpose can be obtained by 
smoothing the data with a median filtering technique. Af- 
ter this operation the histogram is symmetric around its 
mode and the background noise standard deviation can 
be estimated by means of a gaussian fit to the histogram 
itself. 



3.4. Correlation coefficient 

False detections, associated to noise spikes or residual PSF 
features of bright stars, are recognized on the basis of their 
low correlation coefficient with the PSF, which represents 
a template for each true star in the image. The correlation 
coefficient (Gonzalez & Woods 1992) is defined as 



c{a,h) = 



Y.X.V [* (2^' y) ~A[p{x''a,y~b)- p] 



J2x,y [i (^' y) - €JY.x,y [Pi^-a,y~b)- pf 



(2) 



where i{x,y) and p{x,y) are the object and the PSF re- 
spectively, I and p are the corresponding mean values. 
Maximizing the correlation coefficient as a function of the 
offset (a, h) yields an objective measure of similarity. Af- 
ter maximizing c{a^ b) for integral offsets, it is possible to 
repeat the procedure for sub-pixel shifts, improving the 
positioning accuracy. 

The correlation coefficient is computed on the core of 
the PSF: typically the central spike of the diffraction pat- 
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tern is considered, out to the first dark ring. A fair cor- 
relation threshold must be fixed in order to discriminate 
and reject unlikely detections, without losing faint stars 
contaminated by the background noise; a value of 0.7 or 
0.8 is acceptable in most cases. 

The correlation coefficient represents also an effective 
tool to select the stars with the highest photometric re- 
liability, since generally a very high correlation value is 
associated to resolved single sources. 

3.5. Saturated stars 

Saturated stars provide precious information on the PSF 
halo, so it may be useful to include them, appropriately 
reconstructed, in the PSF extraction process. In addition, 
the repaired saturated stars can be recognized by the star 
detection algorithm and their contribution be taken into 
account during the analysis of near and fainter sources. 

The core of a saturated star is replaced with a shifted 
scaled replica of a preliminary estimate of the PSF. The 
repaired star is defined as 



.{x,y) 



f iix,y) 



iii{x,y) < T 



\ fp{x -xo,y- Ho) otherwise 



(3) 



where T is the upper linearity threshold of the detec- 
tor. The position of the center {xo,yo) is estimated by 
means of the maximization of the correlation coefhcient, 
which is not sensitive to the intensity levels. The scaling 
factor / is calculated with a least squares fit to the wings 
of the saturated star. Of course the saturated pixels are 
excluded from all the computations. The background is 
temporarily subtracted in order to prevent affecting the 
repair process. 

This procedure has been applied to the brightest star 
(IRS 7) of the Galactic Center image shown in Sect. 4.2. 
This source is not saturated in the original data, but it 
has been artificially corrupted by an upper cut at half 
maximum. The repair procedure is able to reconstruct the 
original peak with an error of ~5%, imposing a positioning 
accuracy of 1/2 pixel. 



astrometry and photometry are referred to the absolute 
centering and normalization of the PSF array. The opti- 
mization of the parameters is performed by minimizing 
the least squares error between the data and the model. 
If the noise is known, it is possible to weigh the data by 
their inverse standard deviation and obtain a statistically 



optimal fit (Beck & Arnold 1977). In this case, formal 
error estimates on the parameters can be obtained (Bev- 



ington & Robinson 1994). The optimization is performed 



by means of an iterative Newton - Gaus s technique with 
linearized Hessian (Beck & Arnold 1977, Luenberger 1984 
). The M oore- Penrose generaliz ed in verse conce pt (Ru st 
& Burrus|l97|, Bcndinelh et al. |1987| , Lorenzutta |l967D is 
applied to invert the Hessian matrix at every step. The it- 
erations are stopped when the parameters approach a sta- 
ble value. The major difficulty of a Newton-like method is 
represented by the computation of the model derivatives 
with respect to the parameters, some of which (stellar po- 
sitions) yield non-linear dependence. For this purpose the 
Fourier shift theorem is applied: 

p{x-Xn,y-yn) = 



FT-^ 



FT [p{x, y)] e-»2'r(«x„+i,y„)/W 



(5) 



where FT represents the discrete Fourier transform op- 
eration, N is the sub-image size and u, v are spatial fre- 
quencies. The derivative with respect to x„ is then 

dp{x -Xn,y -yn) _ 



dXn 



FT-' 

FT-^ 



FT [p [x, y)] e-^2-(-„+«y„)/^ ( _,^ 



2'KU 1 r , 

-«-^ ) FT[p{x-Xn,y-yn)\ 



(6) 



and requires in practice an interpolation of the PSF to 
compute p{x — x„,y — y„) .In principle this interpolation- 
based method can only be applied to Nyquist-sampled im- 
ages and this is currently the main limit of the algorithm. 
A similar technique has been described by Veran & Rigaut 
(p^98|). 



3. 6. Fitting procedure 

A sub-image centered on the star of interest is extracted 
and approximated with the model 

h(x,y) = so{x,y) + ^ fnp{x - Xn,y - yn) 



+&0 + bix + &22/ 



(4) 



where so{x, y) is the fixed contribution of known stars out- 
side the sub- image support, Ng is the number of point 
sources within the sub-image, a;„, y„, /„ are the position 
and flux of the n — th source, p{x,y) is the PSF and bo, 6i, 
62 are the coefficients of a slanting plane representing the 
local background. It should be stressed that the retrieved 



3. 7. Sampling and interpolation 

A band-limited function in one dimension is Nyquist- 
sampled if the step size fulfills the condition 

Ax < i/. (7) 

where fc is the so-called cut-off frequency of the spectrum. 
In this case it is possible to reconstruct the original con- 
tinuous function from a set of equally-spaced samples by 



means of the so-called sine interpolation (Mariotti 1988) 



For a nearly diffraction-limited PSF, the critical sampling 
condition is generally stated by saying that the its FWHM 
must contain at least two pixels. Many accurate and effi- 
cient interpolation schemes exist to perform the PSF shift 
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required by the fitting procedure. Probably tlie most sim- 
ple is a straightforward application of the Fourier shift the- 
orem, but one may also use bicubic splines or the fast sine 



interpolation algorithm described in Yaroslavksy (1997). 
In StarFinder we have applied the cubic convolution in- 
terpolation method (Park & Schowengerdt 1983), which 



approaches very closely the optimal sine interpolation for 
Nyquist-sampled data; this algorithm is implemented in 
the IDL function INTERPOLATE. 

Even though the interpolation of the PSF array is al- 
lowed only on Nyquist-sampled data, the cubic convolu- 
tion technique seems to produce acceptable results even 
with marginally under-sampled images. Tests performed 
on Airy diffraction patterns, with a sampling step twice 
as large as the critical sampling step size, indicate that the 
interpolation-induced oscillations amount to a few percent 
of the image peak, as opposed to 10-20% of other interpo- 
lation techniques like Fourier shift or bicubic splines. 



4. Applications 

4-.1. Simulated field 

The code was first applied to a simulated image (Fig. 0) 
including 1000 stars, placed randomly in the frame and 
with a given magnitude distribution (Fig. o). Each star 
is a scaled copy of a long exposure high-Strehl PSF, ob- 
tained with the ADONIS AO system at the ESO 3.6 m 
telescope. The image is 368 x 368 pixels large (13" x 13") 
and has a stellar density of 6 stars arcsec"^. Photon, read- 
out noise and a background nebulosity, normalized to the 
same fiux of the stellar sources, were added to the image. 
The faintest stars have a peak signal-to- noise ratio of 5. 

We performed a standard reduction of the artificial 
image using the "default" parameters of the method. 

The PSF was estimated by superposing the images of 
the four brightest stars in the field. The retrieved PSF is 
very similar to the true one (Fig. pi). 

Figure o shows the good agreement between the true 
and the observed luminosity function. The sole statisti- 
cally significant discrepancy is in the bin from magnitude 
7 to 8, where ~25% of the stars were lost; the other small 
differences are due to photometric errors which shift some 
objects to a neighboring magnitude interval. The lost stars 
are ~10% of the total number of sources and are gener- 
ally faint: about 90% have magnitude between 7 and 8, the 
rest is in the bin between magnitude 6 and 7. The only 
bright lost star has magnitude ~4.5 and is the secondary 
component of a very close binary, with a separation of 
just 1/2 pixel. Roughly 70% of the lost stars are located 
at a distance < 1 PSF FWHM from the nearest object, 
~20% are on the first diffraction ring of a brighter source 
and only ^^10% are isolated. It should be stressed however 
that ^15% of the lost stars can be recognized by visual 
inspection as faint objects in the halo of the four brightest 
stars in the field, independently of their separation from 




Fig. 4. Synthetic field image with 1000 sources. The PSF 
Strchl ratio is ^40%. The display stretch is logarithmic. 
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Fig. 5. Axial plot of the true PSF (continuos line) and of 
the retrieved PSF (dashed line). 



the nearest source. It is apparent that the blending effect 
and the contamination by the halo of very bright stars do 
account for the lost stars. Note that the number of false 
detections in this simulated field is negligible (1 case out 
of 1000). 
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Fig. 6. Comparison between the true (continuous line) 
and the estimated luminosity function (dashed line). 

The plots in Figures ^ ^ show the astrometric and 
photometric accuracy of StarFinder. About 80% of the 
detected stars have both astrometric error smaller than 
0.1 PSF FWHM and photometric error smaller than 0.1 
magnitudes. The stars with less accurate astrometry or 
photometry are generally faint: ^75% belong to the mag- 
nitude interval from 7 to 8, ~20% to the bin between mag- 
nitude 6 and 7 and only ~5% are distributed in lower bins 
of the luminosity function. These stars have, in ^^45% of 
the cases, a lost source in their immediate neighborhood 
within the first diffraction ring of the PSF. The others are 
faint stars dispersed in the halo of the brightest sources 
or in the noisy background nebulosity. 

After discussing the performance of the code with a 
standard analysis, it is interesting to examine how the re- 
sults are affected by the main parameters of the method. 
Applying the de-blending strategy we detected ~-|-10% 
more binaries in the separation range between 1/2 and 
1 PSF FWHM, even though the overall detection gain 
is less than 1% referred to the total number of sources. 
With a higher number of iterations of the main loop we de- 
tected ~-|-30% more binaries in the range from 1 to 2 PSF 
FWHM. Decreasing the detection threshold from 3 to 1 
times the noise standard deviation, we found '^+25% more 
binaries in the range between 1/2 and 1 PSF FWHM, but 
with 10 faint (mag > 8) false detections instead of 1. In- 
creasing the threshold on the correlation coefficient, from 
0.7 to 0.8, we reported no false detection, but the number 
of lost stars increased by about 60%; the additional lost 
sources were generally fainter than magnitude 7, but not 
necessarily in crowded groups. Lowering the correlation 



Fig. 7. Plot of astrometric errors vs. relative magnitude 
of detected stars; the errors are quoted in FWHM units 
(1 FWHM ~ 3.6 pixel) and represent the distance be- 
tween the calculated and the true position. A tolerance 
of 1 FWHM has been chosen to find the coincidences be- 
tween the detected stars and their true counterparts. The 
small errors for the brightest stars are due to blending 
effects. 
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Fig. 8. Plot of photometric errors vs. relative magnitude 
of detected stars. The brightest star in the field has mag=0 
by definition. 
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threshold to 0.6 we detected more faint isolated stars and 
binaries, at separations between 1 and 2 PSF FWHM, but 
with a higher probability of false detections (2 instead of 
1). Finally the astrometric and photometric accuracy ap- 
proaches a stable level after a few (^^2) re-fitting itera- 
tions. 



4-2. Galactic Center 

The code was run on a 15 min exposure time K band 
(2.2//m) image of the Galactic Center (Fig. 0), taken with 
the PUEO AO system on the 3.6m CFH telescope (Rigaut 
et al. |1998| ). The Strehl ratio in the image is '--45%. 
0.13", with a samphng of 0.034" 



The PSF FWHM is 
/pixel. The Adaptive Optics gui de sta r was a 77iii=14.5 
star (called star 2 in Biretta et al. 1982 ) located about (to 
the upper left) 20" from the center of the image, out of 
the field of view of the figure. There is therefore a slightly 
elongation of the PSF towards the direction of the guide 
star. However, due to the fact that the isoplanatic patch 
was much larger than the 13" x 13" shown in the figure, 
a space-invarient PSF fits the data very well, as we will 
show. 

A standard analysis was performed, analogous to the 
one described in Sect. 4.1 for the synthetic stellar field. 
About 1000 stars were detected, with a correlation coeffi- 
cient of at least 0.7; the reconstructed image is shown in 
Figure |ol 



We evaluated the accuracy of the algorithm by means 
of an experiment with synthetic stars. We created 10 
frames adding to the original image 10% of synthetic 
stars at random positions for each magnitude bin of the 
estimated luminosity function. The 10 frames were an- 
alyzed separately. As in the simulated case, a distance 
tolerance of 1 PSF FWHM was adopted to find coinci- 
dences between the detected stars and their true coun- 
terparts. The lists of detected artificial stars were merged 
together and the astrometric and photometric errors were 
computed and plotted as a function of the true magnitude 
(Figs. ^ |l^). The plots show no apparent photometric 
bias and high astrometric and photometric accuracy: the 
stars brighter than magnitude 5, for instance, have a mean 
astrometric error < 0.5 mas and a mean absolute photo- 
metric error < 0.01 mag. It should be stressed however 
that the artificial sources are contaminated by the back- 
ground noise present in the observed data and by the pho- 
ton noise due to neighboring stars, but no additional noise 
was added. Figure O shows a comparison between the 
mean luminosity function retrieved in the 10 experiments 
and the truth. Assuming an expected error for each bin 
equal to the square root of the corresponding number of 
counts, according to the Poisson statistic, the only signif- 
icant differences occur for magnitudes fainter than 9. It is 
also interesting to consider the magnitude distribution of 
the false detection cases (dashed-dotted line in Figure |l3|) , 
i.e. the detected stars for which we found no counterpart in 
the original list, within a distance of 1 PSF FWHM. The 
mean percentage of false detections in the 10 experiments 
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Fig. 11. Plot of astrometric errors vs. relative magni- 
tude of detected synthetic stars; the errors are quoted in 
FWHM units (1 FWHM -- 4 pixel) and represent the dis- 
tance between the calculated and the true position. 

is 2% of the total number of stars. The false detections 
are almost always very faint (mag > 8); their number is 
comparable to the square root of the total counts only 
in the last magnitude bin, for magnitudes fainter than 10. 
The percentage of false detections reported in these exper- 
iments seems to confirm the analysis performed by visual 
inspection on the stars detected in the original frame. 

5. The StarFinder code 

The StarFinder code has been provided with a collection 
of auxiliary routines for data visualization and basic image 
processing, in order to allow the user to analyze a stellar 
field, produce an output list of objects and compare differ- 
ent lists, e.g. referred to different observations of the same 
target. The input image is supposed to be just calibrated. 
The code is entirely written in the IDL language and 
has been tested on Windows and Unix platforms support- 
ing IDL v. 5.0 or later. A widget-based graphical user 
interface has been created. The main widget appearing 
on the computer screen is nothing more than an interface 
to call secondary widget-based applications, in order to 
perform various operations on the image. The basic doc- 
umentation about the code can be found in the on-line 
help pages. IDL users might wish to run interactively the 
StarFinder routines, without the widget facilities: com- 
plete documentation on each module is available for this 
purpose. 



Fig. 12. Plot of photometric errors vs. relative magnitude 
of detected synthetic stars. The brightest star in the field 
has mag=0 by definition. 

6. Conclusions and future developments 

The elaboration of real and simulated data seems to prove 
the effectiveness of StarFinder in analyzing crowded iso- 
planatic stellar fields characterized by high Strehl ratio 
PSFs and correct sampling, reaching in this case the full 
utilization of the data information content. The code can 
be applied also to low Strehl data with results compara- 
ble to those attainable by other methods. Moreover we are 
evaluating its performance on undersampled images. 

StarFinder is also reasonably fast: the analysis of the 
Galactic Center image (368 x 368 pixels, ^^1000 stars) re- 
quires between 5 and 10 minutes on a normal PC (Pentium 
Pro - 64Mb RAM - 350MIIz). The graphic interface makes 
it accessible to users unfamiliar with IDL. 

In the near future we plan to provide the code with 
tools to handle images with spatially variable PSF and 
new elements which might be helpful in the analysis of a 
stellar field. An interesting application, still in progress, 
aims at cleaning mixed fields from the contamination of 
foreground stellar images, leaving out true or suspected 
diffuse sources. 

The StarFinder package and its technical documenta- 



tion are directly available at: http://www.bo.astro.it 



or can be obtained on request to E. Diolaiti. 
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